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Abstract. We discuss the signature of a pecuiiar constcliation of very massive 
stars at a projected distance of 2-3 pc around R136a. We discuss various scenar- 
ios for its possible origin, such as independent clusters, triggered star formation, 
supernovae, and ejections via dynamical interactions. If the latter scenario were 
the correct one this would have significant implications on the way to probe 
the conditions in dense stellar cores, and on the evolution of massive clusters in 
general. 



1. Observational findings 
1.1. Introduction 

The 30 Doradus region in the LMC is the largest and most massive H II region 
in the Local Group at a distance of 53 kpc. With in a dia meter of 15' (200 pc) 



it contains more than 8 x 10^ of ionized gas ((Kennicu tt 1984). The inner 
60 pc contain the stellar cluster NGC 2070 with about 2400 OB stars, which are 
respon sible for abou t 1/3 of the total ionizing radiation of the entire 30 Doradus 



region (|Parker 19931 ) . The dense center of NG C 2070 is called R136 (HD 38268). 



With a bolometric luminosity of 7.8 x lO^L© ()Malumuth Sz Heap 1994^ in the 
inner 4.5 pc R136 is often considered the closest example of a starburst region. 

1.2. Observational peculiarity 

Figure^shows the projected distribution of emission line stars in the central part 
of NGC 2070. For the purpose of this study their exact spectral classification 
is rather subordinate; what is important here, however, is that these stars are 
very massive O-stars (see Table HI) . The stars labeled in Fig. ^ either have 
been classified as Wolf-Ravet stars (iParker 19931 ) or X-ray sources, presumable 



HMXBs (|Portegies Zwart et al. 200 
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Figu re 1. HS T image of R136, the central portion of the 30 Doradus region 
(^Massev fc Hunter 1998.1 . The exposure time was 26 s with fiher F336W (cen- 
tered on 3344A a nd with a 38 l l bandpass) using WFPC2. The Wolf-Rayet 
stars identified by ( iParker 19931) are indicated by squares, and the eUipses give 
the locations of the Chandra X-ray sources. The siz es of the ellipses give th e 
la positional accuracy for the Chandra sources l|Portegies Zwart et al. 2002^) . 



Table^lists the projected distance and spectral types for these stars. About 
one third of them is located within 4" of R136a, close to another third is spread 
across the entire region, with distances ranging from 16" . . . 135". But more 
than one third of these massive stars in 30 Dor appear to be located in a shallow 
sphere of r ~ 8" — 13" (although even more of the massive stars, which we 
attribute to the core region, may be actually located in the sphere behind or in 
front of R136a). 

Could this be j ust a statistical fluctuation? For a core radius of Vc ~ 1" 
( Brandl et al. 199 ^1 and under the assumption that the stellar density distribu- 
tion follows roughly a standard cluster King profile, we would expect 
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2.2 stars. 



in the projected sphere, while we find 9 (±3) emission line stars - well above a 
random statistical fluctuation! Thus we conclude that their apparent location 
must have a physical origin. 
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stars within and far from R136 


Stars within the "ring" 


Name 


Sd. tvpe 


Distance 


Name 


Sp. type Distance 


Mk30 


03If*/WN6-A 


16" 


< d< 135" 


Mk33 


03If* 


12.8" 


Mk51 


04If 


16" 


<d< 135" 


TV T 1 o O O 

Mk33S 


03III{f*) 


11.2" 


Mk53a 


WN8 


16" 


<d< 135" 


Mk34 


WN4.5 


10.2" 


R145 


WN6 


16" 


<d< 135" 


Mk35 


03If*/WN-6A 


11.9" 


R140N 


WC5+WN4 


16" 


<d< 135" 


Mk37W 


04If+ 


11.2" 


R140S 


WN4.5 


16" 


<d< 135" 


Mk37 


04If+ 




R144 


WN6-A(B) 


16" 


<d< 135" 


Mk39 


03If*/WN-6A 


12.3" 


R146 


WN4+OB 


16" 


<d< 135" 


Mk42 


03If*/WN-6A 


8.1" 


R136al 


WN4.5 




d < 4" 


R134 


WN6 


11.7" 


R136a2 


WN4 




d < 4" 








R136a3 


WN4.5 




d < 4" 








R136a5 


03If*/WN6-A 




d < 4" 








R136a7 


03III(f*) 




d < 4" 








R136b 






d < 4" 








R136c 






d < 4" 








Table 1. 


Emission line stars in and near R136. 


The spectral type 


is based 



2. Discussion of various scenarios 



2.1. The stars were born at their observed location 



Independent smaller clusters Generally, massive sta rs do not form in isolation 
but in clusters and associations. For a Salpeter IMF (I Salpeter 195fTl ) one would 
expect about 65 stars with m > 4Mq for each star of m > 60 Mq. Are the 
massive stars around R136 surrounded by increased stellar density, relative to 
the overall cluster density profile, or do they appear to be isolated? 

Figure El shows the local luminosity functi ons, derived from W FPC2 data 
in the V (F555W) and I (F814W) fiher bands (|Hunter et al. 19961 V for 16 ran- 
domly chosen test fields and the 9 emission line stars. There is no noticeable 
enhancement in stellar density around the peculiar emission line stars that would 
support the presence of smaller clusters. However, we note the possibility that 
smaller associations could have dissolved within only a couple of million years 
to a level, which is below the detection limit. 



Triggered star formation Sequential star form ation has been fou nd to have 
triggered the L MC associations L H9 and LHIO (" Parker et al. 199^ as well as 
N144 and N158 (jLortet k Testor 1988. ) in the vicinity of the 30 Doradus nebula. 
Several sites of ongoing star formation, such as Hodge 301 and R143, are located 
only 3 arcminutes from R136. However, the peculiar emission line star are 
located very close to R136, in a r egion of ve ry little extinction where the gas has 
been mostly cleared out (jBrandl e t al. 199^ ). No gaseous filamentary structures, 
often associated with triggered star formation are present in the inner parts. 
Most importantly, the emission line stars around R136 are likely of the same 
age than the stars within R136. The core of R136 still appears to be the most 
recent site of massive star formation within 30", making a second generation of 
massive stars in its vicinity rather unlikely. 
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Figure 2. Local luminosity functions for 16 randomly chosen test fields (8 
closer to the core and 8 farther out than the "sphere") and the emission 
line stars, within 0.5pc (2"). The number counts were normalized to the 
overall cluster density profile and the results from the 16 test fields were 
averaged. The same approach was followed to derive the average luminosity 
function of stars around the 8 emission line stars in the sphere. Beyond 
V = 20.5 the data become incomplete in the vicinity of a bright source. The 
data are WFPC2 ob servations in the V (F555W) and / (F814W) filter bands 
ijHunter et al. 1996|) . 



2.2. The stars were born within R136 and moved outwards 

Supernova explosions Runaway OB stars, i.e., apparent field stars with high 
relative velocities of > 40 km/s have been known for a long time — a prime 
exampl e are AE Auriga and /x Columba both moving at 100 km/s away from 
Orion ( Blauuw Sz Morgan 195 J ). Unfortunately, relative velocities are difficult 



to determine for the stars around R136. Although it would only take them 
0.3 Myr to reach the observed location from the cluster core, even at a runaway 
velocity as low as 10 km/s the corresponding proper motion would only be 0."001 
over a 20 year baseline. Radial velocity studies are in progress, but complicated 
by the energetic winds from these stars. 

At any rate, some of the emission line stars are 0-star binaries, which 
cannot be explained by supernova explosions. Most importantly, with an age of 
less than 4 Myr, R136 is unlikely to have yet produced such a large number of 
supernovae, and even more so when taking the kinematical age, i.e., the "travel 
time" from the birthplace to the observed locations, into account. 

Stellar interactions One can esti mate that the time scales for global mass seg- 
regation of R136 are about 10*^ yr ( Brandl et al. 19961 ). However, the relaxation 



time for the high mass (m/^) stars alone is significantly shorter than for the low 

mass (mi) stars by a factor ^1 + ~, where cjl and an are the velocity 

dispersions for the low and high mass stars, respectively. Hence, even in the 
simplified analytic framework, which assumes point-like, identical particles and 
no binaries, the relaxation time for the high mass stars can easily be 1 — 2 orders 
of magnitude shorter. 
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Recently, Portegies Zwart et al. (2004) performed N-body simulations of 
a much more realistic situation: 130,000 stars with a mass distribution given 
by a Salpeter IMF, a cluster half-mass radius of 1.2 pc, and an initial binary 
frequency of 30%^. The initial conditions used in the simulations are similar to 
the conditions in R136. The high central density, existence of binaries, finite 
stellar sizes and realistic mass spectrum allow the cluster to undergo sigificant 
dynamical evolution early on - provided that the core density is sufficiently high. 
Figure 01 shows the evolution of the cluster core radius for different concentration 
parameters. In the most extreme case, core collapse happens at less than 1 Myr. 
At those high stellar densities dynamical processes (binary-binary, or binary- 
single star interactions) can account for numerous ejections of massive stars 
from R136. 




Figure 3. Evolution of the cluster core radius with time for different concen- 
tration parameters ( rdrhm)- The densest system s (large W numbers) show 
dramatic variations l)Portegies Zwart et al. 2004() . 



3. Open questions 

Can massive binaries he ejected and become bright X-ray sources? Whether 
an ejected binary can become a luminous X-ray source depends on its orbital 
separation. For a total mass of R136 of 5-10^ M0, a mean stellar mass of 0.5 M© , 
and a half-mass radius of 1.2 pc a binary of mass rrn, = 40 and a binding 
energy Ebind = 2000A;T will have an orbital separation oi d = 222i?0, which is a 
typical separation for an X-ray binary. 

Why is this signature not observed in NG C 3603 ? Alt hough NGC 3603 has 
been referred to a Galactic clone of R136 ( Moffat 1994^ there are significant 



^The binary frequency of stars, in particular fiigfi mass stars, is a very diffi cult problem to attack , 
both observationally and theoretically. Studies in the Orion cluster by (jPreibisch et al. 1999f) 
using speckle interferometry indicate a large number of 1.5 — 4 companions per primary OB 
star - about 3 — 8 times higher than the frequency of low-mass companions. 
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differences. Most re levant here is that NGC 3603, with an age of about 1 Myr 
(iBrandl et al. 19991 ^. is significantly younger than R136. Given the strong de- 
pendency on central density and quick evolution (see Fig|21) a few hundred thou- 
sand years can make a big difference. 

Why should all ejected stars be observed at the same distance? First, the tidal 
radius of R136 is about 21" (jMevlan 1993l 'l. //the stars are still bound to the 
cluster potential they are in elliptical orbits around R136 and most likely to be 
found at the apocenter. Second, most of the dynamical interactions, which have 
led to an ejection of a massive star, have occurred within a narrow time window, 
namely after the core has reached a critical density. Third, even spherically 
uniform distributions tend to appear ring-like in projection. 
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